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ABSTRACT 

We generate the peculiar velocity field for the 2MASS Redshift Survey (2MRS) catalog 
(|Huchra et al.l 12005) using an orbit-reconstruction algorithm. The reconstructed velocities of indi- 
vidual objects in 2MRS are well-correlated with the peculiar velocities obtained from high-precision 
observed distances within 3,000 km s _1 . We estimate the mean matter density to be il m = 0.31 ±0.05 
by comparing observed to reconstructed velocities in this volume. The reconstructed motion of the 
Local Group in the rest frame established by distances within 3,000 km s _1 agrees with the observed 
motion and is generated by fluctuations within this volume, in agreement with observations. Having 
tested our method against observed distances, we reconstruct the velocity field of 2MRS in successively 
larger radii, to study the problem of convergence towards the CMB dipole. We find that less than 
half of the amplitude of the CMB dipole is generated within a volume enclosing the Hydra-Centaurus- 
Norma supercluster at around 40/i _1 Mpc. Although most of the amplitude of the CMB dipole seems 
to be recovered by 120/i _1 Mpc, the direction does not agree and hence we observe no convergence 
up to this scale. Due to dominant superclusters such as Shapley or Horologium-Reticulum in the 
southern hemisphere at scales above 120/i _1 Mpc, one might need to go well beyond 200/i _1 Mpc to 
fully recover the dipole vector. 

We develop a statistical model which allows us to estimate cosmological parameters from the re- 
constructed growth of convergence of the velocity of the Local Group towards the CMB dipole mo- 
tion. For scales up to 60ft -1 Mpc, assuming a Local Group velocity of 627 km s _1 , we estimate 
Qmh 2 = 0.11 ± 0.06 and <7g = 0.9 ± 0.4, in agreement with WMAP5 measurements at the l-tr level. 
However, for scales up to 100/i _1 Mpc, we obtain il m h 2 = 0.08 ±0.03 and as = 1.0 ±0.4, which agrees 
at the 1 to 2-cr level with WMAP5 results. 

Subject headings: dark matter — mcthods:analytical, numerical and observation 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) dipole in- 
dicates that our Local Group (LG) moves with a veloc- 
ity of 627 ± 22 km s" 1 towards the direction I = 276°, 
b = 30° in galactic coordinates. The first attempt 
at making a comparison between the dipole induced 
by the gravitational influence of structures in our Lo- 
cal UnivBrj3eand_thi3 observed CMB dipole was made 
by lYahil et all (ll980T) using the r evised Shapley- Ames 
catalog and by IDavis fe Huchral (fl98l with the CfA 
catalog. This init ial attempt was followed by a great 
number of works (lYahil et al.l 11986: iLvnden-Bell et al.l 



tion (IKocevski fe Ebelindl2006l. IPl ionis fc Kolokotronis 



19891: IRowan-Robinson et all 119901 : fetrauss et alj 119921 : 
Hudson lll993fT More recently, this comparison has been 
made with the ex tended source catal og of the 2 Micron 
All-Sky Survey bv lMaller et"aLl (|2003l ). These studies all 
agree within 10-30° with each other and with the ob- 
served direction of the CMB dipole. Nowadays, it is con- 
sidered unlikely that the Hydra-Centaurus-Norma super- 
cluster at around 40/i~ 1 Mpc could be solely responsible 
for the dipole. A number of authors suggests that one has 
to go at least as far as the Shapley concentration at about 
150/i _1 Mpc in order to fully recover the dipole mo- 
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119981 iBranchini fe Plionis 
for gravity induced dipole, 

iHoffman etaLl 120011 for POTENT analysis of the ob 
served tidal field). However, due mainly to sparse- 
ness of data at very large distances, there is still no 
consensus on the depth of the convergence for the 
CMB dip ole, as some authors argued for a quicker con- 
vergence (lErdogdu et all l2006al. lErdogdu et all l2006b|. 
Lvnden-Bell et all 119891. lLavavl 119871. lLahav et all 119881 



Strauss et alJ 119921 lYahil et all 11989 and lLahav et al l 
19891 ). This problem is even more complicated by the 



fact that we do not know with a high precision the value 
of the linear bias for each of these surveys. This uncer- 
tainty gives some freedom on the interpretation of the 
apparent con vergence of the velocity o f the Local Group. 
For example, IRowan-Robinson et al.l (|2000l ) argued that 
we reached a co nvergence assum ing a low value for the 
biasing whereas iHudson I (|1994f ) argued against with a 
higher value for the biasing. 

The other approach of studying the convergence of 
the velocity of the Local Group is to use sur v eys o f 
peculiar velocities of galaxies. iWatkins et al.1 (|2009l ) 
showed that all the recent surveys give consistently a 
large bulk flow on a lOOfe " 1 Mpc scale. Only the data of 
lLauer fc Postmanl l|1994f ) still remains inconsistent. This 



seems to indicate that the gravity induced dipole should 
not have converged by lOO/i -1 Mpc. 

The common approach to recovery of the CMB dipole 
has been to use linear theory to reconstruct the veloc- 
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ity field from redshift surveys (see e.g . lErdogdu et al.l 
I2006all2006bl ; iKocevski fc Ebelindl200l . Methods based 
on linear theory can suffer from inadequacy in dealing 
with large fluctuations such as the Hydra-Centaurus- 
Norma supercluster, the Perseus-Pisces supercluster and 
the Virgo cluster. These problems may be enhanced by 
th e coupling with redshif t space distortions. Contrary 
to lErdogdu et al.l (|2006bl ). we also test here our recon- 
structed velocity field against distance measurements. 

Here, we apply the Lagrangian method, Monge- 
Ampere-Kantorovitch (MAK), of peculiar- velocity re- 
construction to the 2-Micron All-Sky Redshift Sur- 
vey (2MRS) catalog (iHuchral l2000t IHuchra et all [20051 ; 
lErdogdu et all l2006al : lErdogdu et all l2006bf) and pro- 
duce a 3-dimcnsional map of the velocity field and study 
the Local Group velo city convergence. Compared to 
lErdogdu et all (j2006bf ). we account better for non-linear 
effects in peculiar velocities, though at the price of hav- 
ing to correct for redshift distortion. The method has 
been adapted to work directly with redshifts and allows 
us to go well beyond the linear theory into the non-linear 
regime. The method has been tested previously against 
simulations and mock catalogs and has been shown to re- 
construct reffaHe_p_e£uliar_velocities on scales above 4-5 
Mpc (see lLavaux et al1l2008l and references therein). In 
addition, using this method we can directly generate the 
3-components of the peculiar velocities and hence over- 
come projection effects. 

The reconstructed velocities of the 2MRS galaxies 
are tested against observed peculiar velocities obtained 
from high-precision dist ance measuremen ts within a 
3,000 km s" 1 (3k) radius (|Tullv et al.ll200l . The recon- 
structed and observed velocities are well-correlated with 
small dispersion. We make an independent estimate of 
fi m using this comparison which is in excel lent agreement 
with the results obtained using WMAP5 (jDunklev et all 
l2009t) . The second test consists of comparing the recon- 
structed velocity of the Local Group in the rest frame of 
the 3,000 km s _1 sample to the observed velocity. Our re- 
sults show that the Local Group motion in the rest frame 
of 3,000 km s" 1 is mostly genera ted within this vo lume, 
as indicated by the observations (jTullv et al.ll2008h . 

Having tested our method, we study the origin of 
the CMB dipole. We reconstruct the velocity field of 
2MRS in successively larger radii in order to determine 
whether there is any trend of convergence toward s the 
CMB dipole. Contrar y to lErdogdu et all (l2006bft but 
more in agreement with iPike fc Hudson! (|2005l ). we show 
that the depth of the convergence in 2MRS lies beyond 
120ft." 1 Mpc. Due to severe incompleteness of the 2MRS 
catalog beyond 120ft" 1 Mpc, we can only put a lower- 
limit on the value of the convergence depth. Our method 
allows us to determine the rate at which the dipole 
is approached: less than 50% of the dipole amplitude 
is achieved at around the Hydra-Centaurus-Norma dis- 
tance and at ~ 120ft" 1 Mpc we recover about 87% of the 
amplitude of the CMB dipole but with no evidence for 
convergence in direction. We then compare the rate of 
convergence that we measure to theoretical predictions 
given by ACDM and the linear theory. It appears that 
our results are in agreement with ACDM when we con- 
sider galaxies at distances lower than 60ft" 1 Mpc. Our 
results deviate la to 2a from expectations of the cos- 
mology of WMAP5 for reconstruction of larger radii. A 



problem relative to the bulk flow on a lOOft -1 Mpc scale 
has also been noted recently by IWatkins et all (|2009f ) 
and on a 300ft " 1 Mp c scale bv iKashlinskv et all (|200l . 
IWatkins et al.l (I2009D show t hat, t hough the bulk flow 
found bv lLauer fc Postman! (|1994f) is not in agreement 
with any other peculiar velocity survey so far, there does 
exist an excessively large bulk flow, whose amplitude and 
direction is compatible with our findings. In this analy- 
sis, the Local Group velocity was taken as a constraint 
for the conditional probability density function. We have 
made a parallel analysis by removing this constraint and 
studying the joint probability distribution function in- 
stead of the conditional probability distribution function. 
We have tested that weighing our likelihood analysis by 
the probability of occurrence of the velocity of the Local 
Group does not change our results. 

This paper is organized as follows. In Section [2l we 
introduce the 2MRS catalog and the 3,000 km s" 1 dis- 
tance catalog and we describe the catalog that we make 
by combining these two. In Section [3j we describe our 
method of peculiar velocity reconstruction. In Section 21 
we constrain the density parameter O m by comparing the 
reconstructed velocity field of 2MRS with the observed 
velocities obtained from measured distances within the 
3,000 km s" 1 volume. In Section O we study the recon- 
structed evolution of the velocity of the Local Group to- 
wards the CMB dipole. We then discuss, in Section I6TTT 
the compatibility of this growth with the cosmological 
parameters as given by WMAP5. We develop a Bayesian 
analysis in Section 16.21 which allows us to make an es- 
timate of cosmological parameters using our measure- 
ment on the velocity of the Local Group in different rest 
frames. In Section [7] we conclude. 

2. THE COMBINED CATALOG: 2MASS REDSHIFT 
CATALOG (2MRS) PLUS CATALOG OF DISTANCES 
WITHIN 3,000 km s" 1 (3K DISTANCE CATALOG) 

The 2MRS (jHuchral l2000l : IHuchra et al1l2005h is the 
most uniformly sampled all-sky redshift survey to date 
and its selection in the near infra-red reduces the impact 
of the zone of avoidance. The catalog is based on galaxy 
selection in the i nfra-red from the Two Micron All-Sky 
Survey (2MASS. ISkrutskie et al.ll20M ). The K=11.25 
mag 2MRS is now complete and contains about 23,200 
galaxies. 2MRS provides good coverage to a mean dis- 
tance of 60ft" 1 Mpc and becomes extremely sparse after- 
ward. 

We fill the zone of avoidance, which mainly extends 
over |6| < 5 — 10°, by mirroring galaxies from the 
unma sked reg ion (first mentione d in Lvnden-Bell et al. 
119891 see also lShava et al.l 119951 and iLavaux et al.ll200S 
for details). The mask has been chosen to cover the 
part of the sky within 10° of the galactic plane for 
galactic longitudes I < 30° and I > 330°, and within 
5° of it for other galactic longitudes. That way, we 
avoid the deeply obscured region behind the galactic 
bulge, at the expense of losing, e.g., the Ophiuchus su- 
percluster. The fake galaxies are represented by red 
points in the top panel of Fig. [T] The kinematic fin- 
gers of god (|Jacksonl Il972f ). associat ed with clusters, 
are co mpressed using the algorithm of IHuchra fc Gellen 
(|1982f l. We used a fiducial velocity V F = 1000 km s"\ 
Vq = 350 km s" 1 , Dq = 0.41 ft" 1 Mpc, which corresponds 
to detecting galaxy count overdensities of 5n/n ~ 80 
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Fig. 1. — Combined 2MRS redshift catalog and 3k catalog of dis- 
tances: Top panel: Large solid green circles are the galaxies for 
which a distance measure is available. The zone of avoidance has 
also been filled and is shown by the red belt running across the 
catalog. Bottom panel: angular positions of a few dominant su- 
perclusters in 2MRS catalog with more distant objects shown by 
redder circles. We represented the position of the Ophiuchus su- 
pcrcluster which has not been included in our galaxy sample. 

(| Crook et all 120071 . I2008T) and provide r obust statisti- 
cal pr operties of groups as explained by iRamella et al.l 
(fl997h . Wc assume a constant M/L ratio as it has 
been shown that the masses of big groups of galax- 
ies o nly moderately depend on luminosities (jLin et alJ 
I2004T ). Moreover, we have checked using the virial theo- 
rem on big groups of galaxies (typically M > 10 12 M Q ) 
that this is the best assumption we can make given the 
present data. Thus, for the calculation of the veloc- 
ity of the Local Group, o ur method looks both l ike a 
flux- weighted computation (|Lvnden-Bell et aT1ll989f ) , be- 
cause we depend on luminos ities of galaxies to infer the 
mass, and number- weighed (jYahil et al .1119911 ). because 
we will be dependent on a correction of redshift distor- 
tions. Incompleteness is taken into acco unt by using 
the Sc hech ter luminosity functi on given in I Crook et al.l 
(f2007h and iHuchra eFaLl (|2005l ) for 2MRS. The specific 
values of the parameters are: $* = 1.06 10 _2 /i 3 Mpc~ 3 , 
ill* = —24.2 + 51og 10 hro- We have verified that, as far 
as galaxies of the 2MRS are concerned, there does not 
seem to be any systematic effect in the relation between 
the K20 magnitude and the intrinsic luminosity. 

We chose to execute the reconstruction itself in the 
CMB rest frame. Though it means that redshift dis- 
tortions may be stronger near the observer, this en- 
sures that better boundary conditions are enforced on 
the outer edge of the reconstruction volume. Using the 



Local Group rest frame for large reconstruction volumes 
may lead to the equivalent of the so-called rocket effect 
(Kaiser 1987), which in the context of MAK still requires 
to be properly modeled. Being in the CMB rest frame, 
this effect should be negligible. 

The second observational component is an extended 
catalog of galaxy distances. Information from four 
techni ques has been i ntegr ated: the Cepheid vari- 
able (iFreedman et all [200lh . Tip o f the Red G i- 
ant Branch (|Karachentsev et all 120031: iLee et al.lll993h 
Surfa c e Brightnes s Flu ctuation ( Tonrv fc Schneider! 
119881: iTonrv et all I2001h. and Luminosi tv-Linewidth 
(|Tullv fc Fisherlll977l : iTullv fc Piercj[2000f) methods. In 
all, there are 1791 galaxies with distance measures within 
3,000 km s" 1 (whence the name "3k catalog"); over 600 
of these are derived by at least one of the first three 'high 
quality' techniques. The 3k catalo g of di stances has been 
described in detail bv lTullv etall (j2008f > . 

The combined 2MRS catalog and 3k distance catalog 
are shown in Figure[TJ In the present study, distances are 
averaged over groups because our method cannot mean- 
ingfully recover the velocities on sub-group scales. The 
present mixed catalog is composed of 24,819 galaxies, 
among which 1,126 has been assigned a distance. Dis- 
tances arc assigned to 109 groups of galaxies out of 695 
groups. There are 617 galaxies with measured distances 
but which have not been grouped. 

3. TECHNIQUE: MAK RECONSTRUCTION OF PECULIAR 
VELOCITY FIELD 

The technique that will be implemented here is 
the Monge-Ampere-Kantorovich (MAK) reconstruction. 
This method provides a recipe for reconstructing galaxy 
orbits that is unique to the degree that orbits can be 
described as following str aight lines under su itable co- 
ordinate transformations (jLavaux et~aT1[200l and ref- 
erences therein). In essence, it tries to find the unique 
displacement field that does not produce shell crossings 
and that map an assumed homogeneous initial density 
field to the presently observed density field. Finding 
this displacement field corresponds to finding the solu- 
tion of a Monge- Ampe re equation, or to so lve a Monge- 
Kantorovitch problem (|Brenier et al.ll2003l ). This prob- 
lem, once discretized, is equivalent to searching for the 
minimum of the discretized action 

N 2 

s = S (%■(»•)- x >) • w 

t=l 

which assigns initial comoving Lagrangian positions 
to final comoving Eulerian positions Xj. We fol lowed a 
different technique than in lLavaux et al.l (|2008| ) to ac- 
count for redshift distortions. Modifying the action using 
Zel'dovich approximation may lead to quite large errors 
for the reconstruction of the velocity of the observer. In a 
volume of radius ^lOh" 1 Mpc, the distortions are dom- 
inant compared to the Hubble flow, and introduce extra 
shell crossings as compared to real space coordinates. We 
use a te chnique inspired from earlier attempt using linear 
theory ijPike fc Hud son 2005]). It consists in running an 
iterative algorithm on the smoothed reconstructed pecu- 
liar velocity field. The steps are as follows: 

(1) Wc let i be the iteration step number. 
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Fig. 2. — 2MRS peculiar velocity field: Top panel shows a thin 
slice of the 2MRS peculiar velocity field adaptively smooth on a 
grid with 128 3 sites sampling a cube with an edge of 240 h~ 1 Mpc. 
The velocity field is then subsampled 4 times before being shown. 
The underlying density field has been computed by putting objects 
at their redshift position. Bottom panel: Peculiar velocities of in- 
dividual 2MRS galaxies in a 40 h Mpc slice, centered on the Su- 
pergalactic plane SGZ=0 km s . 2MRS becomes severely incom- 
plete after I20/i Mpc. Hy-Cn stands for the Hydra-Centaurus 
supercluster, Hr for Hercules, Co for Coma, P-P for Perseus-Pisces. 



(2) We run a reconstruction of the peculiar velocities 
on a catalog Catalog(i) of objects. The coordi- 
nates of these objects are assumed to be a good 
approximation of their real space coordinates. In 
that case we minimize the equation (JTJ) . 

(3) We smooth the reconstructed peculiar velocity in 
redshift space with a Gaussian kernel of radius 



2.5ft." 1 Mpc. This allows us to keep the large scale 
flow intact for the correction while alleviating the 
small scale effects for the correction of the distances 
of the objects of the catalog. 

(4) We move the objects back from their redshift coor- 
dinates from catalog Catalog(O) to their presumed 
new real space coordinates using the velocity field 
obtained at the step (3). This new set of position 
allows us to create the catalog Catalog(i+l) of 
objects. 

(5) We loop to point (1) until convergence. 

We set Catalog(O) to be the 2MRS catalog in redshift 
coordinates. This algorithm in practice converges very 
quickly to the the true peculiar velocity field without 
redshift distortions, typically in 2 or 3 iterations. This 
effectiveness counterbalances the huge cost of running 
several MAK reconstruction. 

4. RESULTS: I. 2MRS VELOCITY FIELD AND 
COMPARISON WITH MEASURED DISTANCES: 
ESTIMATION OF Q M 

In reconstructing the velocity field using MAK 
based on priors, we first set O m = 0.258 and 
bias dKaiser fc Lahavi 119891) to 1 as indicated by 
WMAP5 results (IDunklev et alJ2009l) and us i ng 2dF and 
SDSS /WMAP results (ITeemark et all [2001 : ICole et all 
I2005D . We then do a self-consistency check on the 
presumed cosmological parameters by confronted recon- 
structed and observed Local Group motion. This is the 
approach taken in the next Section (Section [5]). The 
velocity of each galaxy in 2MRS is reconstructed using 
these parameters and checked against observed distances 
within 3,000 km s _1 . This test is detailed in Appendix iDl 

In this Section, we take a different approach and leave 
fl m free and then constrain its value by maximizing the 
correlation between the reconstructed and observed pe- 
culiar velocities. This approach allows us to constrain the 
value of the bias parameter. The 2MRS velocities are re- 
constructed using a uniform grid of size 130 3 sampling a 
cubic volume of 260 3 ft _3 Mpc 3 as shown in Fig. [21 The 
motion of the Local Group is obtained using an inter- 
polation based on the adaptive weighting of the peculiar 
velocities of the objects that lie within 4-5/i _1 Mpc from 
us (method detailed in Appendix [Cj) . We have checked 
that increasing the reconstruction resolution does not sig- 
nificantly change the reconstructed velocities. 

We present, in Fig. 02 the result of the comparison of 
observed peculiar velocity field vs. reconstructed pecu- 
liar velocity field in the volume of radius 3,000 km s" 1 . 
Both fields have been obtained using adaptive smoothed 
interpolation on the the line-of-sight component of the 
velocities (observed or reconstructed) of the objects put 
at their redshift position (Appendix [Q . As we are using 
the redshift coordinates, and not the distance-induced 
coordinates, we should be free of the so-called volume 
Malmquist bias. Moreover, the two fields are enforced to 
be smoothed in exactly the same way at each spatial loca- 
tion, the one needing a larger smoothing scale prevailing 
on the other. Finally, we opted for smoothed interpola- 
tion considering that galaxies are fair tracers of the un- 
derlying continuous velocity field. As it is adaptive, the 
smoothing scale is left free, while in practice it remains 
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in the range [3.4,8.7]/i _1 Mpc. 4 The top panel gives the 
result of the comparison of the raw velocities. However, 
it may happen that the bulk flow of the 3,000 km s _1 vol- 
ume is badly reconstructed. To avoid this problem, we 
subtract the reconstructed (observed respectively) bulk 
velocity of the 3,000 km s _1 region from reconstructed 
(observed respectively) peculiar velocities. The compar- 
ison of the resulting fields is given in the bottom panel 
of Fig. [3l The observed and reconstructed velocities are 
now visually well-correlated. To test the c orrelation, we 
use tw o quantities already introduced in lLavaux et al.l 
(2001). First, we define the correlation coefficient 

(^rec.r^obs.r) 
T = I 9 = ^ ' 

y( € rec,r}(«obs,r) 
With ?) rcc ,r = V Tec ,i - (Wrocr) an( i ^obs.r = W obs ,r - («obs,r), 

u roc ,r the line-of-sight component of the reconstructed ve- 
locity field, u bs,r the line-of-sight component of the ob- 
served velocity field. Second, we define the typical re- 
construction error 

O" 2 = ^((Vrcc.r - W bs,r) 2 ) (3) 

We obtain r = 0.73 and a = 65 km s -1 . These 
values shows that our reconstructed velocity field is 
of good quality. Indeed, the standard deviation a is 
rather small compared to the typical extent of the PDF 
(~400km s^ 1 ). The correlation r is better than what has 
been obtained for mock catalogs for which observational 
effects has been included (r mock ~ 0.53), in spite that the 
comparison is done on a smaller number of tracer of the 
peculiar velocity field: 576 tracers within 3,000 km s _1 in 
the case of this paper against ~1600 in mock catalogs). 
This seems to indicate that our resu lts should be as goo d 
as the one obtained on simulations (|Lavaux eFaI][200l . 

From the comparison between reconstructed to ob- 
served peculiar velocities, it is possible to measure fi m . 
Though a pure likelihood ap proach would seem the natu- 
ral way, lLavaux et al.l (|2008D showed the high sensitivity 
to the prior of a likelihood which does not take into ac- 
count correlations of the velocity field. For this paper, 
we will th us use the apparent l y less sensitive met hod of 
moments (jColombi et al.1120071; lLavaux et aLll2008f) to de- 
termine f2 m on the adaptively smoothed velocity field. 
Smoothing has the advantage of using the correlation 
of the velocity field to increase the signal-to-noisc. The 
noise coming from both errors on observed distances and 
from MAK modeling will be in particular diminished. As 
it has been shown in this paper that the slope between 
reconstructed displacements and observed peculiar veloc- 
ities 

/ ( g rec,r) , A s 

s= \jwu (4) 

seems not to be statistically biased, we will use it here 
as our estimator of f2 m . We improve the technique by 
computing a set of different slopes from data points that 
lies farther and farther from the perfect correlation (i.e. 
the diagonal in Fig. 0). This allows us to make an esti- 
mate of the error bar on our result on f2 m through the 

4 This range corresponds to the selected smoothing radius for 
95% of the volume. 



use of the two o ther estimator of th e slope r x s and s/r 
(Appendix B of lLavaux et~ani2008D . 

We assumed h — 0.80, a value compatible with 
the catalog of distances with the current calibration 
of distance indicators. The value assumed for h only 
m arginally affect s our e stimation of f2 m as it was shown 
in lLavaux et al.l ([2008). From the 576 velocities that 
were reconstructed, we esti mate Q m = 0.31 ± 0.05 . This 
result is in agreement with lPike fc Hudso n ( 2005). This 
is only consistent at the lower end of f2 m with previ- 
ous results based on orb it reconstruction methods (e.g 
IMohavaee fc Tullvll2005f ). An extensive likelihood anal- 
ysis including correlations, which would yield a better 
estimation, is postponed to another paper. 

We estimate a systematic error of ~ 9% on recon- 
structed peculiar velocities due to the assumed values 
of cosmological parameters. To this error, we add in 
quadrature a random reconstruction error of 70 km s _1 
according to the mean (both on amplitude and by compo- 
nent), as estimated using reconstruction on simulations 
in redshift space (Fig. Appendix [B| . 

One must be aware of the d anger of cosm i c vari ance, 
as it has been highlighted in lLavaux et alJ (|2008| ). A 
way of checking it is to use the density field of the whole 
2MRS as an indicator of the local density fluctuation. We 
computed a density contrast pzo/pwo — 1 = —0.046 for 
the Universe within 30/i -1 Mpc, compared to the mean 
density within lOO/i" 1 Mpc using the mass density field 
obtained from 2MRS data. It means the volume on which 
we make the comparison is slightly undcrdense and thus 
may bias our value to lower r2 m . However, the amount 
of the systematic bias due to cosmic variance seems to 
be here negligible compared to errors due solely to the 
noise. Nevertheless, this value is in good agreement with 
other measurements such as the one given by WMAP 
(|Dunklev et all 120091 ) for the ACDM model. Having 
obtained and tested our reconstructed velocities within 
3,000 km s _1 , we go to larger scales and study the origin 
of the motion of the Local Group with respect to CMB, 
Vlg/cmbj m the next section. 

5. RESULTS: II. 2MRS VELOCITY FIELD AND THE 
ORIGIN OF THE CMB DIPOLE 

The CMB dipole motion, is obtained by using the re- 
constructed 3d velocity field generated using the 2MASS 
Redshift Survey. The velocities are reconstructed within 
increasing radii centered on the Local Group. A table of 
the resulting dipole as a function of the reconstruction 
radius is given in Table [1] The entries arc represented in 
Fig. 0J No convergence is achieved below 120/i _1 Mpc. 
The 2MRS sample becomes highly incomplete beyond 
this scale and further conclusions cannot be made before 
a more complete sample becomes available. Whether the 
Shaplcy supercluster yields the convergence to the CMB 
dipole also remains questionable, as farther away in the 
southern hemisphere, dominant structures such as the 
Horologium-Reticulum supercluster could change the di- 
rection of the dipole. Thus, it is possible that the depth 
of the convergence towards CMB dipole lies well-beyond 
the Shapley concentration itself. 

We also plot the increase in the bulk flow relative to 
the CMB dipole over the reconstructed peculiar veloci- 
ties within 120/1" 1 Mpc in the bottom panel of Fig. [2j 
The figure demonstrates that less than half of the am- 
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Fig. 3. — Observed vs. Reconstructed velocity field — These two 
panels give the result of the comparison between the reconstructed 
and the observed line-of-sight component of the velocity field. The 
top (respectively bottom) panel indicates the result of a compari- 
son before (respectively after) having corrected for bulk flows. We 
have used the method of moments using the part of the scatter 
within the 95% probability contour in both cases. For the top 
panel we obtain O m ~ 0.47. For the bottom panel, we estimate 
f2m ~ 0.31. The correlation coefficient r in that case is equal to 
0.73, with a typical a posteriori velocity error of 65 km s" 1 . To 
plot these two distributions, we assumed h = 0.80. Though the 
centroid of the scatter has been repositioned for the estimation of 
the slope, this correction has not been applied for the above distri- 
butions. Decreasing h would move the scatter to the right of the 
diagonal. The filled contours limit the regions that corresponds to 
50% (black), 68% (dark gray), 95% (gray) and 99% (light gray) of 
joint probability to have a reconstructed and an observed velocity 
at some position in the Local Universe. 



plitude of the dipole is generated within the volume en- 
closing Hydra-Centaurus-Norma. To reach convergence 
a significant contribution to the dipole has to be made 
from the Shapley supcrcluster and beyond. This agrees 
with an analysis of the X-ray data (jKocevski fc Ebelingl 
I2006T ) . We also fail to observe convergence by the time 
we reach 120/i _1 Mpc. Indeed, although the amplitude 
of the reconstructed dipole seems to approach that of 
CMB, its direction remains well away from it. Using ob- 
servational error bars, we estimate that the observational 
error on the direction of the Local Group velocity should 
not exceed ~ 6° at 95% confidence. In addition, we esti- 
mated from an iV-body ACDM simulation that there is 
an intrinsic 22° error (95% confidence) in the direction 
of the reconstructed velocity because of the modeling er- 
rors of the MAK reconstruction. In sum, we expect the 
observed and reconstructed vectors to agree to within 
~ 23° whereas they lie at ~ 40 — 50° from each other. 
This problem is highlighted by Fig. where we have 
represented the angular separation between the recon- 
structed Local Group velocity and the CMB dipole. We 
note that beyond ~ 60/i _1 Mpc the separation remains 
at roughly 40 — 50°. We now proceed to compare our 
results to predictions given by the ACDM model. 

6. COMPATIBILITY WITH THE ACDM MODEL 

The measurements of the convergence of the velocity 
of the Local Group to the CMB dipole can be used to 
check the ACDM model. Indeed, the statistic of peaks 
and voids in the density field imprints the fluctuation 
and the rate at which the velocity of the Local Group 
converges to its expected final value, which should be 
the one given by the CMB according to the current cos- 
mological paradigm. We will use here the cosmolog- 
ical p arameters estimated by WMAP5 (IDunklev et all 
20091) and weak lensing measurements (|Fu et al.l 120081 ; 
Benjamin et al.ll2007t iMassev et al.ll2007f) . We use a sta- 
tistical mo deling of the Local Group velocity growth de- 
veloped by [Juszkiewicz et al.1 JT990) (equation 12) and 
ILahav eta l. (1990). The prediction is based on linear 
theory and is able to give the probability of occurrence 
of our meas urements given the cq smol ogical model. A s 
proposed by IDunklev et all (|2009D and lFu et all (|2008D . 
we use fi m = 0.256 (the density of baryon is fib = 0.044 
and the density of cold dark matter is f2 c = 0.212) 
and erg = 0.77. We also set the Hubble constant to 
H = 71.9 km s _1 Mpc -1 and the spectral index n = 1. 
In addition to the cosmology specified above, we enforce 
that the amplitude of the velocity of the Local Group is 
627 km s _1 . 

This conditioning allows us to avoid the problem of the 
probability of our own velocity in the different universes 
we consider in this work. Thus, though we may be at 
a special point in the Universe, our results should not 
be affected by the fact that the velocity of Local Group 
may be unusual. In the following, we thus focus on the 
conditional probability of obtaining a set of reconstructed 
velocities of the local group {v^} given a specification 
of the velocity of the Local Group Vlg- Nevertheless, 
we have also checked in Section HT21 that using the joint 
probability of having both {v^} and Vlg only slightly 
changes the probability contours in the specific case of 
our velocity as given in Table [1] This indicates that 
the velocity of the Local Group is not too uncommon in 
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Fig. 4. — Velocity of the Local Group in progressively larger rest 
frames ,Vlg : The observed amplitude and direction of the CMB 
dipole motion are shown by the horizontal cyan band in the lower 
panel and the solid cyan large square in the top panel, respectively. 
The lower panel shows the amplitude of the velocity of the Local 
Group in successively larger rest frames as the velocity field of 
2MRS is reconstructed at increasing radii. Incompleteness is illus- 
trated by the black error bar on data points beyond \20h~ 1 Mpc in 
the bottom panel. The top panel shows how the direction changes 
as the radius increases. The red curves in the bottom panel indi- 
cate the prediction of growth of the velocity of the Local Group 
for a WMAP5 cosmology. The solid curve gives the expectation of 
the reconstructed velocity for a survey whose radius is indicated 
by the X axis. The two dashed curve indicates the lc fluctuation 
relative to the expectation given by the model. To compute these 
curves, we used the WMAP5 parameters: the density of cold dark 
matter Q c = 0.2 12, the density of baryo ns O b = 0.044, h = 0.719, 
erg = 0.77 and a Eiscnstcin & Hu (1998) power spectrum (without 
baryonic wiggles). 

universes whose cosmological parameters are selected by 
the likelihood analysis. 

We use a ACDM power spectrum as given by 
lEisenstein fc Hul (|1998D but without incorporating bary- 
onic wiggles. We have checked that introducing wig- 
gles does not change the prediction much though the 
introduction of baryons does decrease the expectation 
of the reconstructed Local Group velocity for distances 
<60/i _1 Mpc. This behavior is expected as baryons tend 
to suppress density fluctuations below the sound hori- 
zon while they are linked to photons by the mean of 
the Compton effect. The size of the horizon at the mo- 
ment when ba ryons separate from p hotons is typically 
~ 45/i _1 Mpc (|Eisenstein fc Hulfl99cl . which is the same 
scale at which we observe a change due to the introduc- 
tion of baryons. In the presence of baryons the density 
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Fig. 5. — Misalignment angle of the reconstructed Local Group 
velocity — We have plotted here the misalignment angle between the 
reconstructed Local Group velocity and the direction indicated by 
the CMB dipole. The symbols and colors that are used here are 
the same as in Fig. [4] The red curve represents the 95% probability 
limit of the misalignment for a ACDM universe whose parameters 
have been chosen as estimated by WMAP5. The horizontal black 
thick line gives the expected misalignment, at 95% of probability, 
between the reconstructed velocity and the observed motion of the 
Local Group. It has been estimated by applying the reconstruction 
to one ACDM simulation. 

field has less power on smaller scales, so it is more dif- 
ficult for the Local Group to acquire its velocity using 
only small scale fluctuations. Fully computing the ex- 
pected value of the Local Group velocity for a given sur- 
vey depth, we indeed observe that it decreases by 5-15% 
when we take into account baryons in the power spec- 
trum. 

Now, we may also consider the effect of changing a%. 
Its principal effect is to change the amount of fluctua- 
tion of the velocity field around its expected value. A 
growth of convergence that is slow and regular corre- 
sponds most likely to a low local as, whereas a growth 
with a lot of independent fluctuations favors a high erg. 
Its impact on the expectation of the amplitude of the 
velocity field is more complicated. Indeed, cosmologies 
with high as tends to have stronger fluctuations relative 
to the expected velocity, which yields an higher expected 
amplitude. Thus higher as should increase slightly the 
expectation of the amplitude. This means that even if 
we use only the evolution of the amplitude of the Lo- 
cal Group we should be sensitive to as, though more 
marginally than on the shape of the power spectrum 
f2 m /i. The impact of as is however dominant concern- 
ing the fluctuations of the direction of the velocity of the 
Local Group. Only universes with a high as allow this 
direction to depart significantly from the one given by 
the CMB dipole at scales larger than GOh^ 1 Mpc as we 
will see in the section UTTl 

f2 m represents the true dynamical content of the uni- 
verse. For a given realization of density fluctuations, the 
dynamics is faster for a high r2 m than for a low f2 m . Thus 
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the convergence of the velocity of the Local Group should 
be quicker in a universe with a high fi m than in another 
one with a low r2 m . But the expectation of this conver- 
gence does not tell anything about the probability on the 
velocity of the Local Group itself. We have decided in 
this work not to presume anything on the probability of 
occurrence of the velocity of the Local Group but only 
to check the consistency of the dynamics of the Local 
Universe. 

We are now going to use the prediction of the growth 
of the Local Group velocity as both a checking of the 
ACDM model and a way to estimate cosmological pa- 
rameters. We will explore a limited number of param- 
eters: the total mean density of matter f2 m , the mean 
baryon density f2b, the Hubble constant H and the vari- 
ance us of density fluctuations within a sphere of radius 
8I1- 1 Mpc. 

In Section 16. 1[ we compare our results to predic- 
tions using WMAP5 cosmological parameters. In Sec- 
tion !6.2[ we estimate cosmological parameters using both 
our measurements and the constraint coming from nucle- 
osynthesis. 

6.1. Comparison to WMAP5 results 

First, we assume the ACDM model described above 
and we obtain the red curves represented in the bottom 
panel of Fig. [4j The solid curve gives the expectation 
of the convergence of the dipole given our motion. The 
dashed curves give the estimated la fluctuation for the 
expectation of the velocity of the Local Group for a given 
survey depth. This expectation is computed assuming 
that the Local Group is moving at 627 km s _1 rela- 
tive to the CMB, as determined using CMB experiments. 
Though it seems astonishing to be systematically below 
the solid curve such a behavior may be explained by the 
fact that we are dealing with correlated measurements 
(see Section f6 . 2 1) . For each point, the ACDM model 
seems to be in agreement with the Local Group velocity 
reconstructed by MAK at the la level. A similar predic- 
tion can be obtained for the maximum angle separation 
between the CMB dipole and the reconstructed velocity 
(equation 13). This prediction is represented by the solid 
red and thick curve in Fig. For a given survey depth, 
the misalignment has a 95% probability to be located 
below the curve. This is in perfect agreement with the 
results given by the MAK reconstruction until a depth of 
60/i~ 1 Mpc. There, the misalignment goes up to ~ 50° 
and stays there. At lOO/i -1 Mpc, we have a probability 
less than 3 x 10~ 5 that the misalignment is greater than 
50°. 

Given the importance of the disagreement between the 
predicted maximum misalignment and the observed one, 
we see three options. First, the reconstructed velocity 
of the Local Group may be exceptionally plagued by 
rcdshift distortions, even though on the rest of the vol- 
ume the reconstructed velocity field is very well recon- 
structed. In this case, we are in the sa me position as 
number - weighed methods, which like in lErdogdu et alJ 
(|2006ah yield a large misalignment. Doing a complete 
Monte-Carlo analysis on mock catalogs of these prob- 
lems are beyond the scope of this paper, and we post- 
pone this to future work. Second, we may miss some rel- 
atively important structures in the masked region. Wc 
remind that due to severe incompleteness wc arc miss- 



ing information on structures below 5° for |/| > 30° 
and 10° for |/| < 30°. One of the brightest X-ray clus- 
ters of the Local Universe is the Ophiuchus Cluster at 
8400 km s^ 1 and it may well reside in a hidd en super- 
cluster of significance (Wa kamatsu et al~ll2005f ) . Second, 
more dramatically, it would mean that large bulk flows 
exist at a scale of 100ft. -1 Mpc, typically ~ 500 km s _1 . 
This is extremely unlikely in a ACDM cosmology as 
it is already highlighted by Fig . \5j A similar result 
has been obtained recen tly by iFeldman et alJ (|2008l ) 
and iWatkins et alJ (|2009f ) but using observed velocity 
fields. On the contrast, using linear theor y on 2MRS 
(lErdogdu et al.ll2006at lErdogdu et al.ll2006bl) gives a low 
value of b ulk flow and d isagre es with IFeldman et ahl 
dSH and IWatkins "eTaLl p009Tl . 

6.2. Parameter estimation 

From the measurements of the growth of the Local 
Group velocity, it is possible to estimate the joint proba- 
bility of Q m and as by a maximum likelihood approach. 
To do that, we use a likelihood analysis on the velocity of 
the Local Group in the different rest frames that we ob- 
tained using 2 MRS in Fig.lH A simila r analysis was first 
introduced by iKaiser fc Lahavl (|1989f ) to compare pre- 
dicted dipoles to observed bulk flows. Here, we focus on 
the expected evolution of the bulk flow on different scales 
and how it should compare to our own reconstructed 
bulk flow in a given cosmology. We detail, in Sec- 
tion l6.2.ll the statistics of the growth of the Local Group 
velocity using linear theory , an id ea originally intro - 
duc ed bv Uuszkiewicz efaLT (119901), lLahav et all (|1990f ) 
and lStrauss et alT i 1992f )~In Section [6.2.2[ we see how 
to mix observational errors into this statistical modeling. 
Finally, in Section I6.2.3[ we see the results of this analy- 
sis on the whole set of data points that we obtained on 
the value of the cosmological parameters. 

6.2.1. Statistics of the growth of the Local Group velocity 

To estimate the likelihood of the cosmological pa- 
rameter vector p, we need the conditional probability 
-P(vi, V2, . . . , vat|Vlg7P) of reconstructing the velocities 
of the Local Group {v^} in the set of rest frames {IZi} 
given that the true velocity of the Local Group is Vlg 
and the cosmological parameters are described by p. To 
estima te this probability, we will follow lJuszkiewicz et al.1 
(|1990f ) and use linear theory and assume the components 
of the velocity field are independent Gaussian random 
variables. To ease the notation, we write ilk to designate 
the N dimensional vector whose i-th component is the 
fc-th component of the i-th vector: itfe,, = (vj)fe. Us- 
ing linear theory, we can compute the variance Ci of the 
velocity of the Local Group computed using a homoge- 
neously sampled spherical survey of depth Ri 
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P{k)W 2 (kRi)dk 



k=0 



with W{x) 
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sin(x)/x. We also define 
{(3Hf 



6tt 2 



P(k)dk 



(5) 



(6) 



fc=0 



the variance of the velocity field. Similarly, one may 
obtain the correlation coefficient between the com- 
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poncnt k of Vj and Vj 
1 



■{Vk,iVk,j) 



P{k)W{kR % )W(kR 3 )dk (7) 

&-K z Uiaj J k=0 

We also need the correlation I\ between Vf~,i and Vfc.LG 
1 



6ir 2 <7iav Jk= 



P{k)W{kRi)dk. (8) 



fe=0 

We note that, even in the non-linear regime, the veloc- 
ity field remains highly Gaussian. However, these co- 
efficients do not take into account potential distortions 
to its correlations due to the MAK reconstruction and 
non-linear effects. Though, we expect the computation 
of these coefficients to be correct on large scales, pre- 
cise measurements using TV-body simulations would be 
required here to improve the precision on the final es- 
timation of cosmological parameters when smaller scales 
are used. This measurement could be interesting as small 
scales may contain more information on Qb and ag. Us- 
ing these quantities it is now possible to compute the 

— > 

probability P(vi, v 2 , . . . , v N \ Vlg) = -P(^|V L g), with 



it = (Ui, il2,it3). We refer the reader to the appendix 
for the details of the co mput ation. We carry here the 
result given in equation (|A10[) : 



exp f-i Y. W^i&JkjM.j^j 



with 



ilk,i FiVk y LG 

JjJk,i = — ■ 

Oi ay 



(9) 



(10) 



and M s the top-left part of the inver t of the covariancc 
matrix as defined in appendix by Eq. (jAlj) . As the three 
components of the Local Group velocity arc independent, 

the total probability of measuring the tensor H given 
that the velocity of the Local Group is Vlg is thus: 



p(Ii|v lg ) = f[P(il k \V kthG ) . 



(11) 



k=l 



However, in Fig. [3J we see that errors coming from re- 
constructions are quite significant. We thus have to take 
them properly into account into the likelihood analysis. 



J 



6.2.2. Measurement errors 

To account for reconstruction errors, we assume that they are independent from one survey depth to another and 
they are properly modeled by a Gaussian distribution. We will thus write the probability that the component k of the 
reconstructed velocity v rcCi j of the Local Group for a survey cut i, given that the true velocity should be for some 
cosmological model and the expected error dispersion for the component k is a e> k,i, 



-P^reCjfcjil^fcjij fe,fc,z) 



(2^/v fc / xp U 



(12) 



where u rcc ,fc,i = (v rec ,i)fc and Vk : i = (vj)^. To ease the notation, we use now the iV-dimensional vector 2Jfc, whose 
components are QJj jt = (v reCj j)fe. So these vectors holds the component k of the reconstructed velocities for all survey 
depth. 

We may now merge the two probabilities (J9j) and (fT2|) to obtain the total probability to measure QJfc given a cosmology 
determined by the vector of parameters p and the component k of the velocity of the Local Group Vlg,A; : 



P(%k\V hG ,k,p) = j^dU (j[P(v hG ,i,k\ U i^e,i,k)) P(Ui,U 2 ,...,U N \V hG , k ,p), 



(13) 



with U a vector with N components. One may use a similar transformation as used in (|A4[) to compute the above 
integral and obtain the sought above probability: 



P(Kk\VL Gik ,& e ,k,p) 



(2^W detAW6 e , fc ) ^ VhG ^ &e ' k) 



(14) 
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where we defined the following quantities 



N / s 9 N 

/ $n. _ tt 



m v,e)=j2( g mcaiM ) + E m™ a n,iW~,Ms^ , (15) 

i=l ^ 1 ' i,j=l 

2Ua,ca,u(», V, 6) = CUa»,i(9J, V, 6)/^ - T^/dy , (16) 

j=i V j ay fc=i / 



2 c 

MobM^SHM.Ai + ^y*. (18) 

It is interesting to note that f/ m ean,j is a vector with A~ components which corresponds to the best average velocity 
given the cosmological model and observations. In some sense, it is a velocity whose observational noise has been 
subtracted thanks to the modeling. Finally, we may assemble the statistics of the three independent variable 53^ with 
k = 1,2,3. We may do so by multiplying their respective probability distribution given by Eq. (|14[) to obtain 

N 

P^, . . . ,v N \V hG ,(T e ,i, . . . ,a e , N ,p) = YlP(fOk\VuG,k,e k ,p). (19) 

fc=i 

We now use this probability and the Bayes theorem to compute the probability for the cosmological parameters p 

P(P Vi,...,V W ,V L G,ffe,lr--,« r e,JV = r , „, pTJ 7 (20) 

J p dp P(Vi , . . . , Vjv | V lg , CT e ,l , • ■ • , (Te,N,P)P(P) 

with P(p) a prior on the probability of p. In practice, we use a flat prior on £l m h 2 and os parameters and impose a 
prior on ilbh 2 and P as we detail in the next section. 



6.2.3. Results 

We set here a Gaussian prior on the Hubble constant 
to h = 0.72 ± 0.08 as specified by the HST Key project 
(jFreedman et al.ll2001h and the spectral index to ns = 1. 
Nonetheless, we checked that the result depends only 
weakly on this prior. We also fix the observed veloc- 
ity of the Local Group to Vlg/cmb = 627 km s -1 . As 
we ar e using the power spectrum of lEisenstein fc Hul 
(|1998f ) without wiggles the cosmological parameter vec- 
tor is thus p = (S! m /i 2 , ilbh 2 , as, h). Note that the like- 
lihood analysis, depending on the density power spec- 
trum, should naturally depend on the shape Q m h and 
not on fl m h 2 . On the other hand, WMAP5 data yield 
fi m /i 2 , so we arbitrarily chose WMAP5 parametrization 
over the natural one. The expected growth of the Lo- 
cal Group velocity presents a substantial degeneracy in 
the (tt m h 2 , Qbh 2 ) plane for scales smaller than the sound 
horizon (~ 45/i _1 Mpc). To get better constraints, we 
need to introduce some prior on one of these two variable. 
Followi ng the results obtai ned using the deuterium abun- 
dance (jPettmi et al.ll2008f ). we set a prior on fl^ti 2 to a 
Gaussian distribution centered on u>t, = il^h 2 = 0.0213 
and whose width is 0.0010: 

. , / 1 (uj b - 0.0213\ 2 \ 

p W Kop (-2(-iMirJ J' (21) 

with a constant of proportionality independent of p. We 
are not using the results from WMAP5 to make an in- 
dependent comparison of their constraints and ours. We 
assume a flat prior on ui m = fl m h 2 and as though we 



enforce them to be within the ranges ui m £ [0, 1.0] and 
a 8 G [0.05,7]. 

The results given by Eq. (|20|) are illustrated in 
Fig. El In the left panel, the 68% (respectively 95% and 
99.7%) isocontour of the marginalized joint probability 
of (f2 m /i 2 ,cr 8 ) is represented in filled black (respectively 
dark gray and light gray) contour when all velocities up 
to 100/t. -1 Mpc arc taken into account. In the right panel, 
we represent the isocontours of probability when one con- 
siders only velocities up to 60ft -1 Mpc. The 68% (re- 
spectively 95%) probability for cosmological parameters 
using WMAP5 data alone are given by the hashed dark 
red (respectively hashed light red) filled contours in both 
panels. 

It can be seen that our results in the right panel are 
compatible within la with WMAP5, as far as uj m and as 
are concerned. However, in the left panel, the gray con- 
tours are significantly farther from WMAP5, though still 
within a 2a limit. This behavior is expected as we have 
already noticed that the misalignment angle does not 
seem to be compatible with a standard WMAP5 cosmol- 
ogy, though the growth of the amplitude is compatible. 
The introduction of points beyond ~60/i _1 Mpc shifts 
the contours of probability to low u m and high as- A 
high CTg allows strong fluctuation in the direction, thus 
increasing the expected amplitude, whereas a low f2 m 
limits the impact of these fluctuations on the amplitude 
of the velocity by making the convergence slower. At 
the same time these contours become smaller as the con- 
straints are enforced by a larger number of data points. 
The analysis of all data points up to 100ft -1 Mpc shows 
that such a misalignment is less unlikely than what we 
would have expected from the observation of the curve of 
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Fig. 6. — Marginalized distribution of the joint probability of Q m h 2 /o-g - In red, we represented the isocontour at 68% (dark red) and 
95% (light red) of marginalized probability as given by WMAP5 (CMB data only). Left panel (a): We represented the isocontour at 68% 
(dark gray) and 95% of probability (light gray) as given by the likelihood analysis on the growth of the velocity of the Local Group up to 
lOOfe -1 Mpc. For th e gray contours, we introduced a prior on the mean density of baryons obtained by measuring the deuterium abundance 
UPettini et al.ll2008h . which gives U b h 2 = 0.0213 ± 0.0010. Right panel (b): The dark (68%) and light (95%) gray contours 

have been obtained in the same way as the gray contours of panel (a) except that we limited the depth of surveys to 60/i Mpc inclusive. 



Fig. [5l This can be explained by the fact the misalign- 
ment angles are not independent random variables. A 
more independent set of random variables would be the 
difference between the acceleration of the velocity the 
Local Group when one goes from one survey to a deeper 
survey. Thus, the strong misalignments of ~ 50° for all 
R > 60/i -1 Mpc, as we have on Fig. might be gener- 
ated by a single rare event within 60 — 100/i _1 Mpc from 
the observer which is not included in our present cata- 
log, say a major supercluster in the Ophiuchus direction. 
The relative absence of matter in our catalog compared 
to the reality would push away the reconstructed veloc- 
ity from the CMB direction, which corresponds to what 
we observe. We need to observe in greater detail the 
obscured part of the sky in these range of distances to 
resolve this issue. Another related reason for this mis- 
alignment might be an excessive mass in our model in 
the Perseus-Pisces supercluster. In the region delimited 
by a sphere of 15/i _1 Mpc centered on z=5000 km s -1 , 
I = 140.2, b = -22, we have p ut 7.45 x 10 15 M Q . Ac- 
cording to IHanski et al.l (|2001D , this is on the high end 
of the expected mass of the Perseus-Pisces but not too 
excessive. A lower mass for this supercluster would push 
the velocity of the Local Group in the direction of the 
CMB. 

We may check the adequacy of the conditional prob- 
ability P(23|Vlg, of reconstructing the set of Lo- 
cal Group velocities 2J given our observed motion Vlg 
with the motion of the Local Group itself. To do that, 

we look at the total probability P(2J, Vlg| S ,p) = 



P(23|V LG , e,p) x P(V LG |p). We checked that in this 
case our results are not significantly affected, slightly 
moving the contours of both panels of Fig. [5] towards 
higher u> m and higher as- This behavior is expected. This 
probability takes into account the fact that it is more 
likely to observe the velocity of the Local Group with re- 
spect to CMB in a denser universe and with larger den- 
sity fluctuations. However, these considerations are in- 
sufficient to significantly change the computed contours. 

On the contrary, we may relax the condition on the 
asymptotic reconstructed velocity, which is in principle 
given by the observed CMB dipolc. In that case, we 
rely only on the evolution of the reconstructed velocities 
with the depth of the survey. Such an analysis shows the 
growth of the reconstructed velocity of the Local Group is 
compatible within la. However, this agreement is mainly 
due to the extension of the contours of probabilities of 
Fig. [6] because we do not enforce the observed motion of 
the Local Group. Doing so, we becomes less sensitive to 
potential large scale bulk flows and are unable to draw 
any meaningful conclusion about the cosmology. 

To summarize, our results of this section indicate that, 
on one hand, the total matter density might be lower in 
the local Universe than that measured using WMAP5. 
On the other hand, the amount of fluctuation, erg, is 
more likely to be higher. Quantitatively, the maximum 
of the likelihood for the gray contours of the right panel 
is located at uj m ~ 0.08 and erg ~ 0.57. The mean of the 
gray distribution of the right panel are given by f2 m /i 2 = 
0.11 ±0.06 (Q m ~ 0.21 ±0.11) and cr 8 = 0.90 ±0.42. For 
the gray contours of the left panel, this position is shifted 
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to ft m h 2 = 0.08 ± 0.03 (O ro - 0.15 ± 0.06) and a s = 
1.0±0.4. The most likely parameters are u> m — 0.090 and 
<t% = 0.58 for the right panel, and w m = 0.070 and erg = 
0.78 for the left panel. We note that putting our estimate 
of f2 m of Section U and the above estimate of uj m , we may 
compute ft ~ 0.51, which is at odd with the estimate from 
the HST Key project. We expect some ov erlap at the low 
oj m an d high as part of the contour with Wa t kins et al.l 
((2009h . These discrepancies may be explained in three 
different manners. First, some more work on simulations 
are required to understand the potential systematics on 
these particular measures. Second, we know that there 
exists a large concentration of galaxies in the direction of 
the O phiuchus at about 80ft -1 Mpc (jWa kamats u et all 
|2005[) . This concentration is estimated to be as large as 
the Hercules supercluster and would need to be properly 
accounted in our statistical modeling as it is likely to 
systematically affect our reconstructed velocity. Third, 
we may live in a very special part of the Universe and 
we would thus need to go to larger scales using deeper 
catalogs to clarify this issue. 

7. CONCLUSION 

We have evaluated the 3-dimensional 2MRS peculiar 
velocity field using a Lagrangian method (MAK) which 
works well into the nonlinear regime, on scales above 
4-5ft -1 Mpc. The method has been adapted to work di- 
rectly with redshifts. The reconstructed velocities are 
wcll-corrclatcd with the distances observed within a ra- 
dius of 3000 km s -1 (3k sample). The reconstructed 
velocity of the Local Group in the rest frame of the 3k 
sample also agrees well with the observed value. In ad- 
dition to this test, we have successfully compared the 
observed velocity field to the reconstructed ones within 
this same volume. This comparison has lead to an esti- 
mate of O m = 0.31 ± 0.05. 

We have then studied the origin of the Local Group mo- 
tion in the CMB rest frame by using our 3d reconstructed 
velocities. We have shown that less than half of the CMB 
dipole could be generated within a radius enclosing the 
Hydra-Ccntaurus-Norma supercluster. We have demon- 
strated how the trend of the convergence varies up to 
120ft -1 Mpc and have shown that convergence in po- 
sition to the CMB dipole is not reached even by this 
distance. 

We checked that our measurements of the amplitude 
arc qualitatively in good agreement with the ACDM 
model for cosmological parameters given by WMAP5. 
We also note a weak dependence of the theoretical ex- 
pectation of the amplitude on the presence of baryons 
that could be of interest to get better constraints if fu- 
ture work manages to reconstruct sufficiently well (at less 
than ~5% error) the local dynamics. However, misalign- 
ment angles are significantly larger than anticipated by 
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the ACDM model on scales larger than 50-60ft -1 Mpc. 
To quantitatively check our measurements, we developed 
a Bayesian analysis of the growth of the velocity of the 
Local Group in the context of the linear theory but in- 
cluding correlations. We have shown that the growth of 
the velocity of the Local Group may be a powerful inde- 
pendent tool to explore cosmological parameters and of 
the dynamics of the Local Universe. This analysis yields 
an independent measurement of (f2 m ft 2 , a^) assuming the 
quantity of baryons using the nucleosynthesis theory. For 
scales up to 60ft- -1 Mpc, our measurements agree with 
WMAP5 at the la level. For scales up to 100ft -1 Mpc, 
this agreement drops to 1 to 2a. This problem may have 
different origins: the incompleteness of the 2MRS, bad 
correction of redshift distortions at the introduction of 
the Perseus-Pisces or a real large bulk flow on the scale 
of ~100ft -1 Mpc, which can represent a challenge for 
ACDM . More work on mock catalogs are needed to clar- 
ify this issue, which is beyond the scope of this paper. 

Forthcoming deeper and more complete redshift sur- 
veys especially those in X-ray, could soon establish 
whether one has to go well beyond 120ft -1 Mpc to re- 
cover the CMB dipole or whether the Shapley concentra- 
tion at around 150ft -1 Mpc is sufficient to finally reach 
the convergence. Such a larger survey would also allow 
a better comparison of our analysis of observations with 
the ones made on simulations on a similar scale to check 
A CDM features in the large scale velocity field, such as 
in lStrauss et all (|1995f ). 
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TABLE 1 

The reconstructed CMB dipole in progressively larger rest frames R 



Observed velocity: V LG /cmb = 627±22 km s _1 , I = 276 ± 3, b = 30 ± 3 
"LG/CMB,! = 56 ± 28 km s _1 , u l g/CMB,j = -540 ± 25 km s -1 , ti LG/CMBiS = 313 ± 30 km s~ 



R 






Vlg/r 












v x 


Vy 




1 V| 


I 


b 


angular 


{h' 1 Mpc) 


(km s — x ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(deg) 


(deg) 


separation 


20 


-157 ± 71 


-83 ± 70 


73 ± 70 


192 


207 


22 


62 


30 


-122 ±71 


-173 ±72 


166 ± 72 


269 


234 


38 


35 


10 


-137 ± 71 


-285 ± 74 


196 ± 72 


372 


244 


31 


27 


50 


-172 ±72 


-270 ± 74 


177 ± 72 


365 


237 


29 


34 


60 


-236 ± 73 


-282 ± 74 


101 ± 70 


382 


230 


15 


45 


80 


-208 ± 72 


-230 ± 73 


182 ± 72 


360 


228 


30 


41 


100 


-328 ± 76 


-277 ± 74 


197 ± 72 


473 


220 


25 


49 


120 


-378 ± 78 


-319 ±76 


228 ± 73 


544 


220 


25 


19 


150 


-257 ± 74 


-413 ± 79 


133 ± 71 


504 


238 


15 


38 



Note - Glossary of the symbols used in the above tables. I and b are the galactic longitude and latitude respectively. R is gives the radius 
of the sub-volume of the catalog on which the reconstruction is achieved. V{ are the Cartesian galactic coordinates of the reconstructed 
velocity of the Local Group. The x axis points towards (I = 0,6 = 0), the y axis to (I = 90°, b = 0) and z to b = 90°. 
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APPENDIX 
APPENDIX A 



CONDITIONAL MULTIVARIATE GAUSSIAN RANDOM VARIABLES 

Wc consider N + 1 Gaussian random variables described by the vector (U, V) with U a vector of dimension N and 

V a scalar. In this appendix, we compute the statistics of U given V. We put the normalized variables U{ = Ut/at 

and V = V/ay- We assume that the correlation coefficients 7^ = (UiUj) between the components Ui and Uj of U 

are given. We also assume that I\ = (UiV) is the correlation coefficient between V and {/,. The covariance matrix of 
this TV + 1 variables may thus be written as 

'1 if i = j; 
7i,j ifi<N and j < N ; 



Ti ifi<Nandj = N + l; ^ 



Tj if t = N+ 1 and j < N. 

We now write the invert of this matrix M = C _1 and M s the top-left most N x N sub-matrix of M. In other words, 
this sub-matrix corresponds to the first N lines and TV columns of M. The probability of the vector (U, V) is thus: 

P(U, V) = (2 ^ )(jV+1 ) /2VH ^ *P (-^W-M,,^^, - \m (n+iUN+1) V^ . (A2) 

Now we use Bayes theorem to compute P(U\V): 

P W V ' = /TO7' (A3) 

Before computing the integral in the denominator, we need to rewrite the argument of the exponential: 

*UM S U + 2VM it(N+1) Ui + M (N+1)AN+1) V 2 = 'WM.W + aV 2 . (A4) 

with 

W = U + A, (A5) 

Ai = (M^) h3 M N+ht V = -C N+1 ,iV = -T t V , (A6) 
a = M N+ i. N+1 - M N+ i^M s .ijM N+lj . (A7) 

Wc may now compute the denominator: 



We may rewrite the argument of the exponential in the numerator using (|A4j) and deduce the conditional probability 

P (U|V) = exp (—(Ui - r 4 V)(^ - I .1 ;.W_ ; ) . (A9) 

This conditional probability is written in terms of the normalized variables U and V. The probability may be rewritten 
in terms of the original random variables by an appropriate change of variables: 

N 



F ( U l y ) = (ft v) CXP (-5^/^ - riV/wm/ai - TjVMM.,^ . (A10) 



APPENDIX B 



COMPARISON OF LINEAR THEORY AND MAK RECONSTRUCTION OF DIPOLE 

We have demonstra ted, in contrast to previous reconstruction of 2MRS peculiar velocities (jErdogdu et al.ll2006al ; 
lErdoedu et"al"1l2006bD. th at our method agrees well with recent works which make direct use of observed peculiar 
velocities (jFeldman et all 120081 ). In this appendix, the difference between linear theory and MAK determination of 
bulk flow is studied. 

There are mainly two ways of implementing linear theor y for getting the veloc ity of the Local Group. One method 
is the flux- weighed method (as recently done for 2MRS in lErdogdu et al.ll2006af ). The second method uses linearized 
fluid equation to infer the velocity field from the density field. The advantage of the former method is that it bypasses 
the redshift space distortions, which could be problematic for objects highly extended in redshift space (e.g. Virgo 
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TABLE 2 

Structures probed for comparing with linear theory 



Structure 


^min 

(km s- 1 ) 


^max 

(km s- 1 ) 


I 

(deg) 


b 

(deg) 


Opening angle 
(deg) 


Virgo 


no minimum 


3000 


279 


74 


20 


Hydra 


2000 


4300 


269 


26 


15 


Centaurus 


2000 


4300 


302 


21 


15 


Perseus-Pisces 


4000 


6000 


150 


-13 


20 


Fornax 


no minimum 


1600 


240 


-50 


20 



TABLE 3 

Impact of structures on the final velocity of the Local Group, smoothed with a Gaussian kernel of 5/i -1 Mpc radius 



<5v s ,, 



Method 


Structure (S) 


Mass 


\SV\ 


I 


b 






(M Q ) 


(km s ) 


(deg) 


(deg) 


MAK 


Virgo 


1.3 x 10 15 


89 


302 


88 


Linear 






230 


297 


68 


MAK 


Fornax 


3.6 x 10 14 


52 


261 


-15 


Linear 






94 


271 


-1 


MAK 


Hydra-Centaurus 


4.3 x 10 lb 


88 


281 


35 


Linear 






189 


297 


39 


MAK 


Perseus-Pisces 


8.9 x 10 lb 


105 


167 


10 


Linear 






52 


139 


48 



Note - We are giving in this table the direction and amplitude of the velocity vector corresponding to the velocity of the Local Group 
when all structures are present minus its velocity when we removed the indicated structure. The structures are denned in Table [2] The 
comparison is done for both the result given by MAK and the linear theory. 

cluster). The disadvantage of this method is that it is limited to the velocity of the Local Group. It does not yield 
the peculiar velocity field and especially in the local neighborhood. In addition, to study the impact of farther and 
farther structures on the Local Group, one often substitutes distances by redshifts. 

The second way of obtaining peculiar velocity field from linear theory relics on linearized fluid equations. The 
disadvantage of this method is that the Local Group velocity becomes sensitive to redshift space distortions. The 
advantage is that we obtain the 3d velocity field everywhere and not just at our location. In this appendix, we 
compare this way of implementing linear theory with MAK. 

To estimate the velocity of the Local Group with linear theory, we follow the same procedure for accounting for 
redshift distortions as for MAK, except that we use the linear relation between peculiar velocity field and the density 
field given by: 

divv = -/3H6, (Bl) 



where = fl 



m 9 /b, b is the bias, H the Hubble constant and 5 the density fluctuations. Wc use a Gaussian smoothing 
filter of 2.5/i _1 Mpc of radius and increase progressively (3 to increase numerical stability (e.g lPike fc Huds"onll2005f) . 

To study the difference between linear theory and MAK, we look at four most prominent structures (listed in 
Table [2]) in 2MRS catalog and study their effect on the Local Group motion. We evaluate the Local Group velocity in 
the presence and in the absence of these four structures using both linear theory and MAK. We computed the velocity 
of the Local Group based on a 60/i _1 Mpc cut of the 2MASS Redshift Survey and in the rest frame of this volume. 

Table [3] demonstrates the effect of these four structures on the Local Group velocity. We note that linear theory 
overestimates the Local Group motion due to lower mass nearby objects and underestimates the effect of high mass 
farther objects. This deserves a full theoretical understanding in the framework of perturbative theory which is beyond 
the scope of this work. Since Perseus-Pisces is the most massive object in our study, we also show in Figure [7] the 
effect of this structure on the entire velocity field, as reconstructed using linear theory and MAK. The large scale flows 
are similar in the two methods and, as expected, the small scale flows are different. 

We have done a more extensive test of the difference between linear theory and MAK reconstruction through a Monte- 
Carlo test. We have taken fifty observers, randomly located in the simulation, and built a mock catalo g centered on 
each o f these observers. These mock catalogs are extracted from the ACDM simulation described in lLavaux et alj 
(|2008D . Each mock catalog is limited in depth to 60/i _1 Mpc and mock galaxies and groups are put in redshift space. 
We calculated the velocity of the mock observer given by the linear theory and by MAK and then compared these to 
the true value given by the simulation. The result is shown in Figure [8] MAK presents a dispersion ^30% less than 
linear theory. (For other comparisons with linear theory see Lavaux 2009.) 
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Fig. 7. — Reconstructed velocity field in the region of the Perseus-Pisces in Galactic coordinates - We give here the reconstructed velocity 
field using cither linear theory (left panels) or by solving the Mongc-Ampere-Kantorovitch problem (right panels) for two cases: with the 
Perseus-Pisces supercluster (top row), and removing the galaxies of the Perseus-Pisces supercluster (bottom row). In all cases, we limited 
the reconstruction to a volume of radius 60/i~ 1 Mpc from the observers. The slice is centered on the Perseus cluster and does not contain 
Milky Way. The slice corresponds actually to GY=25/i -1 Mpc. 
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Fig. 8. — Comparing MAK to linear theory through Monte Carlo tests: reconstructed velocity vs. simulated velocity - Wc represent 
here the comparisons between reconstructed velocities (top row MAK and bottom row linear theory) and simulated velocities of 57 mock- 
observers in 57 mock catalogs. The velocity fields (simulated, reconstructed by MAK, reconstructed by linear theory) were smoothed to 
2.5h~ 1 Mpc with a Gaussian kernel. Left panels: direction misalignment vs amplitude of the simulated velocity. Right panels: amplitude 
of the reconstructed velocity vs amplitude of the simulated velocity. The mock-observer velocities reconstructed using MAK (top row) 
show ~30% less dispersion with respect to linear theory. 
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APPENDIX C 



ADAPTIVE SMOOTHED INTERPOLATION 

We make use of a technique inspired from adaptive smoothing (jColombi et al.|[2007f ). while changing the weight 
such that we make an smoothed interpolation of the required field (in this work the velocity field). We put A(y) the 
underlying continuous field, and At its value at the position x^. We now define the smoothed interpolated field 



to-m^iwrh (C1) 



with M the number of tracers on which the smoothed interpolation is done, 

M 

N(y) = J2W(y-Xi) (C2) 

i=l 

the normalization coefficient of the smoothing and 

!1 - \x 2 + jX 3 X < 1 
1(2 -xf Kx<2 (C3) 

otherwise 

the smoothing kernel used. In practice, in the above smoothing procedure, we set M = 32, and R(y) is half the 
distance to the 32-th neighbors from y. 

It can be shown that for class 1 function 5 and an unbiased distribution of tracers of A(y) this smoothing procedure 
converges normally on any volume V. Locally, the sum (|C1[) may be expanded in 

, at / v _ \ 

A W = ^ £ W w W + ( v ^)(y) • - y)) ( C4 ) 



N(y) f-[ \ R(y) 

( M i 

:A(y) + (VA)(y)-j^2 M 



%)nw i(x -- y)| - (c5) 

The second term on the right-hand-side may be bound by: 

( V ^)(y)' (E^^(l^)( x >-y)) <2\\VA\\ +00 , V MR. (C6) 

As R cx n -1 / 3 with n the number density of tracers and that ||V-i4||+oo is bounded on any volume V, the above 
quantity decreases to zero when n — > +oo and the interpolation A converges to A. 

APPENDIX D 



TEST OF 2MRS RECONSTRUCTED VELOCITIES WITH OBSERVED DISTANCES WITHIN 3,000km S" 1 

The reconstructed velocities of objects lyi ng only within the 3,000 km s _1 radius can be compared to the measured 
distances given by the 3k distance catalog (jTullv et al.N2008l ). The measured distances give a velocity of the Local 
Group with respect to the 3,000 km s _1 volume of Vlg/31< = 302±22 km s _1 , I = 241 ±7, b = 37±6. The observation 
indicate s that most of thi s velocity comes from the push from the Local Void and the gravitational pull of the Virgo 
cluster (jTullv et alJl2008T ). The velocity of the Local Group with respect to the 3,000 km s _1 volume (VLG/3k) is 
obtained using our reconstructed 2MRS velocities. The reconstructed V^q/^ is compared with the observed value in 
Fig. [5J The coordinates of the reconstructed dipoles are given in Table 0J The influence of other structures, like the 
Hydra-Ccntaurus-Norma at 40-50/i _1 Mpc, seems marginal in this rest frame. So we conclude that the reconstruction 
indeed shows that the Vlq/3^ motion seems mainly generated within the 3,000 km s _1 volume. 



Class 1 functions arc continuous and their derivative is also continuous. 
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Fig. 9. — The origin of motion of the Local Group in the rest frame of galaxies in the 3,000 km s (3k) volume, Vlg/3]c : The observed 
amplitude and direction of the V l g /3k are shown by the horizontal cyan band in the lower panel and the solid cyan large square in the 
top panel, respectively. The observational error on the direction is shown on the top panel by the thick cyan ellipse. Each reconstructed 
point is affected by a 22° error. The reconstructions indicate that most of the motion is generated by the structures within this volume, in 
agreement with the observations. It is seen that most of the motion is generated within 3,000 km s~ l . The jump at around \20h~ 1 Mpc 
is most probably caused by the severe incompleteness of 2MRS at this distance. The represented error bar have been computed by adding 
in quadrature the 9% systematic error bar due to the uncertainty on C m and a random reconstruction error of 70 km s -1 . Incompleteness 
results in uncertainties illustrated by the black error bars on data points beyond 120/t" 1 Mpc in the bottom panel. 



TABLE 4 

Comparison between the observed and the reconstructed velocity of the Local Group in the 3,000 km s _1 rest frame 





Observed velocity: 


^LG/3k = 302 km s 


~\ « = 241, fe = 37 












VLG/3k,x = -125 ± 24 km s \ « LG /3k,y 


= -210 ± 21 kin s" 


> ""LG/Sk.z 


= 181 ± 18 km s- 1 










^LG/3k 












Vx 


Vy 




\v\ 


I 


b 


angular 


(h.- 1 Mpc) 


(km s ) 


(km s ) 


(km s ) 


(km s" 1 ) 


(deg) 


(deg) 


separation 


30 


-122 ±71 


-173 ±71 


166 ± 71 


269 


234 


38 


5 


10 


-150 ±71 


-228 ± 73 


166 ± 71 


320 


236 


31 


7 


50 


-183 ±72 


-187 ±72 


145 ± 71 


300 


225 


28 


16 


60 


-232 ± 73 


-182 ±72 


90 ± 70 


309 


218 


17 


28 


80 


-183 ±72 


-128 ±71 


154 ± 71 


271 


215 


34 


21 


100 


-282 ± 73 


-176 ±72 


152 ± 71 


371 


213 


24 


27 


120 


-311 ± 73 


-230 ± 73 


164 ± 71 


121 


216 


23 


25 


150 


-231 ± 73 


-234 ± 73 


124 ± 71 
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225 
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